Anomalous Doping Effects on Charge Transport in Graphene Nanoribbons 
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We present first-principles calculations of quantum transport in chemically doped graphene 
nanoribbons with a width of up to 4 nm. The presence of boron and nitrogen impurities is shown 
to yield resonant backscattering, whose features are strongly dependent on the symmetry and the 
width of the ribbon, as well as the position of the dopants. Full suppression of backscattering is ob- 
tained on the 7r — 7r* plateau when the impurity preserves the mirror symmetry of armchair ribbons. 
Further, an unusual acceptor-donor transition is observed in zig-zag ribbons. These unconventional 
doping effects could be used to design novel types of switching devices. 
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The ability to single out a single graphene plane, 
through an exfoliation process [l| , or by means of epitax- 
ial growth Q , has opened novel opportunities for explor- 
ing low dimensional transport in a material with remark- 
able electronic properties [3[. Additionally, the devel- 
opment of graphene-based nanoelectronics has attracted 
much attention owing, to the promising large scale inte- 
gration expectations @, 0] • However, with 2D graphene 
being a zero-gap semiconductor, its use in an active elec- 
tronic device such as a field effect transistor (FET) re- 
quires a reduction of its lateral size to benefit from quan- 
tum confinement effects. Graphene nanoribbons (GNRs) 
are strips of graphene whose electronicproperties depend 
on their edge symmetry and width |5|], and can be ei- 
ther patterned by plasma etching @, 0] or derived chemi- 
cally Q . Band-gap engineering of GNRs has been exper- 
imentally demonstrated Q, and GNRs-based FET with 
a width of several tens of nanometers down to 2 nm have 
been characterized 

Chemical doping aims at producing p-doped or n- 
doped transistors, which are crucial for building logic 
functions and complex circuits [|[. Doping also allows 
new applications such as chemical sensors, or electro- 
chemical switches [Io|. In carbon-based materials p- 
type ( ?> type) doping can be achieved by boron (nitro- 
gen) atom substitution within the carbon matrix [ill ]. 
For metallic carbon nanotubes (CNTs), Choi and co- 
workers reported that boron (B) and nitrogen (N) impu- 
rities yield quasibound states which strongly backscatter 
propagating charge for specific resonance energies [l2j |. 
The interplay between those resonance energies and ex- 
ternal parameters (electric or magnetic fields) may also 
enable the design of novel kinds of CNT-based switching 
devices [13|, whereas the existence of quasibound states 
related to topological defects was unveiled by STM mea- 
surements [14j . 

In this Letter, we report on an ab initio study of the 
effect of both p-type (B) and n-type (N) doping on quan- 



tum transport in GNRs with widths within the experi- 
mental scope. In contrast with CNTs, doping in GNRs 
turns out to display more complex features depending on 
the dopant position, ribbon width and symmetry. Theo- 
retically, two types of GNRs with highly symmetric edges 
have been described, namely zig-zag (zGNRs) and arm- 
chair (aGNRs) [5]. Some recent works have reported on 
the effect of doping in extremely narrow zGNRs [Hj], 
mostly with a width in order of 1 nm. 

Following Ref. [16], we refer to a zGNR (aGNR) with 
N zig-zag chains (dimers) contained in its unit cell as a N- 
zGNR (N-aGNR). We have studied aGNRs and zGNRs 
with widths between 2.3 and 4.2 nm, already within the 
current experimental scope [1, [l7| . The scattering poten- 
tial around the dopant is obtained using first principles 
calculations (SIESTA code [18[) within the local density 
approximation 

We start with the case of aGNRs. Ab initio studies 
show that aGNRs are always semiconducting Ji| 23. 24] 
with width-dependent bandgap scaling. We have studied 
the 20-, 34- and 35-aGNRs, with widths of 2.3, 4.0 and 
4.2 nm, respectively. In carbon nanotubes, two accep- 
tor (donor) quasibound states have been predicted below 
(above) the charge neutrality point (CNP) at low energy 
values when a single carbon atom is substituted by a 
boron (nitrogen) impurity 12, 25[. 

However, in contrast with CNTs, the energies of the 
quasibound states in GNRs are strongly dependent on 
the position of the impurity with respect to the GNR 
edges. A clear increase in binding energy of the bound 
state associated with the broad drop in the conductance 
in Fig. [I] is observed as the dopant approaches the edge. 
The large variation of resonant energies with dopant posi- 
tion indicates that random distribution of impurities will 
lead to a rather uniform reduction of conductance over 
the occupied-states part of the first conduction plateau. 
This is in sharp contrast to the case of CNTs, where 
resonant energies do not depend on the position of the 
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FIG. 1: (color online). Left: Conductance of the 35-aGNR as 
a function of energy for different B positions. The dashed lines 
correspond to conductance for the undoped case. Insets: same 
as in main frame for two selected nitrogen dopant positions 
(at the edge (E), and off-center (B)) for the 20-aGNR. Right: 
unit cell of the 35-aGNR showing the considered dopant po- 
sitions. The carbon atoms are shown in grey and the passi- 
vating H atoms in pink; the other colored atoms represent B 
or N atoms, each color referring to a corresponding colored 
solid line conductance curve. 



dopant around the tube circumference. Our results for 
the 34-aGNR (not shown) confirm this behavior also for 
semiconducting aGNRs. In addition, for certain dopant 
positions, symmetry effects yield a full suppression of 
backscattering even in the presence of bound states. In- 
deed, when B is placed at the exact center of the ribbon, 
the transmission at the first plateau is found to be insen- 
sitive to the presence of the impurity (Fig. [I]- left, bottom 
curve). Conversely, as the defect approaches the ribbon 
edge, the energy resonance floats up towards the CNP 
and the conductance is clearly degraded by the quasi- 
bound states induced by the impurity (Fig. [I]- left). 

To understand such phenomena, we must first note 
that GNRs, unlike CNTs, do not always present a well 
defined parity associated to mirror reflections with re- 
spect to their axis. An ideal odd-index aGNR retains a 
single mirror symmetry plane (perpendicular to the plane 
of the ribbon and containing the ribbon axis), and its 
eigenstates will thus present a well-defined parity with 
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FIG. 2: (color online). Conductance as a function of the 
energy (left) and bandstructure (right) of the 35-aGNR [2(| . 
Dashed lines correspond to the undoped ribbon; solid lines 
correspond to the case of B at the center. The solid blue 
(dashed red) arrow shows the first (second) band below the 
charge neutrality point (CNP) for the doped case. 



respect to this symmetry plane (see Figs. [3}} and [3^). 
The eigenstates of the doped ribbon keep the same par- 
ity with respect to this mirror plane provided that the 
potential induced by the dopant preserves this symme- 
try [l2|, [2?| • For the case of an odd-index aGNR, this 
can only occur when the dopant is located exactly at the 
central dimer line, as is illustrated in Figs. [3fci and [3]i. 
Comparing the wavefunction at the T point associated to 
the first band below the CNP for both the ideal ribbon 
(Fig. [3t>) and the doped (B at center) one (Fig. [3^), 
we observe that states around the CNP are only weakly 
affected by the impurity, indicating that there is no mix- 
ing with neighboring bands (which present an opposite 
parity). As a result, backscattering does not occur de- 
spite the presence of impurity states with energy values 
within the first plateau (Fig. O dashed red arrow, and 
Fig. [3]i). For any other position of the dopant, the well- 
defined parity of the wavefunctions will not be preserved, 
as is shown in Figs. [3b and [3f for the B off-center case 
(red atom (B) in Fig. [1]). In this case, coupling between 
all states restores backscattering efficiency, yielding a full 
suppression of the single available conduction channel at 
a certain resonance energy determined by the precise lo- 
cation of the dopant (red curve (B) in Fig. [1]). 

In comparison with B-doping, the impact of N impu- 
rities on the ribbon transport properties manifests in a 
close symmetric fashion with respect to the CNP. This 
is illustrated in the insets of Fig. [TJ where the conduc- 
tances for a 20-aGNR with N off-center (B) and at edge 
(E) are shown. The same symmetry considerations afore- 
mentioned also apply in the case of odd-index N-doped 
aGNRs. 

We consider now the case of doped zGNRs, and present 
two different systems: namely a 12-zGNR (« 2.4 nm 
width) and a 20-zGNR (« 4.1 nm width). zGNRs are 



FIG. 3: (color online), a-c) 2- dimensional projection of the wavefunctions at the F point associated with the first occupied 
band below the CNP (solid blue arrow in Fig. [2]), for a 35-aGNR with B at center (a), no dopant (b), and B-off center(c). d-f) 
Same as in a-c) for the second occupied band below the CNP (dashed red arrow in Fig. [2}. The position of the B impurity 
is marked in green. Blue or yellow color corresponds to opposite sign of the wavefunction. Dashed lines show the ribbon axis 
(located exactly in the middle of the ribbon). 




FIG. 4: (color online). Conductance as a function of the energy for the 12-zGNR (left) and 20-zGNR (right) with a single B 
dopant placed at the ribbon edge (E, top panel) and off-center (B, bottom panel). The dashed lines correspond to the undoped 
case. The red or light blue colored atoms in the unit cell show the position of the dopant for each ribbon. 



known to display very peculiar electronic properties, with 
wavefunctions sharply localized along the GNRs edges 
at low energies, which significantly affect their transport 
properties [28[. Simple tight-binding models found that 
zGNRs are always metallic with the presence of sharply 
localized edge states in the vicinity of the Fermi level [5[ , 
whereas spin-dependent ab initio calculations report on 
a small bandgap opening up 0. Here spin is neglected, 
and we focus on the effect of a boron defect on the trans- 
port properties. 

Fig. 2] presents the conductance for B-doped 12-zGNR 
(left) and 20-zGNR (right). In contrast to aGNRs, where 
the acceptor character of boron is maintained regardless 
of the position of the impurity, B-doped zGNRs exhibit 
both acceptor and, more unexpectedly, donor character 
when the dopant is placed either at the center or at the 
edges, respectively. This effect may be related to the 
competition between two different phenomena, namely 
the Coulomb interaction of charge carriers with the ion 
impurity and correlation between charges at the edges. 
Our results for the narrower 12-zGNR B-doped ribbon 
confirm such prior prediction [30] for N-doped zGNRs 
[3lj . However, for the wider 20-zGNR, a defect located 
at the edge has almost no impact on the conduction effi- 
ciency in the plateau around the CNP. This is in striking 



contrast to the results for the aGNRs, where doping has 
its maximum effect when the impurity is placed at the 
edge. This suggests that modification of the electronic 
properties of zGNRs solely by means of edge doping 
or functionalization might not be significant on zGNRs 
wider than a few nanometers. On the other hand, the 
geometry of symmetric (even index) zGNRs does not al- 
low a symmetry axis going through a single impurity [32[ . 
As a result, the suppression of backscattering observed 
in the case of aGNRs cannot take place here. 

In conclusion, ab initio calculations of charge trans- 
port in boron and nitrogen doped GNRs with a width 
of up to 4.2 nm have been reported. Doping effects de- 
pending on the ribbon symmetry and width were un- 
veiled, such as a full suppression of backscattering for 
symmetry-preserving impurity potentials in armchair rib- 
bons, and upward energy shift for the quasibound state 
resonances in both armchair and zigzag ribbons. Such 
predictions could be experimentally confirmed, follow- 
ing recent progress achieved on the STM exploration 
of graphene edg es 1331 ] and on boron-doped oriented py- 
rolytic graphite [34|. Finally, chemical doping could be 
used to enlarge the bandgap of a fixed GNR width, re- 
sulting in the enhancement of device performances [351 ] . 
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